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SUmXELY 

The   olsject  of  this  experiment  was  to  investig?.te  the 
extent  of  the  mixing  region  of  a  pulsating  air  jet  at  varying 
frequencies  in  a  steady  secondary  air  flow.  This  was  accom- 
plished by  maicing  a  velocity  survey  of  the  region.   The  study 
v;as  made  \inder  isothermal  conditions  at  three  different  fre- 
ouencies. 

The  results  of  this  investigation  show  that  the 
pulsating  flow  mixes  more  rapidly  and  has  a  shorter  potential 
core  than  the  steady  flov;  jet  of  the  same  configuration. 
Mixing  is  the  best  and  the  potential  core  is  shortest  for  the 
pulsating  flow  at  resonant  frequency.  Decay  of  the  pulsation 
amplitudes  along  the  centerline  is  linear.   Decay  of  the  pul- 
sation amplitudes  laterally  is  non-linear.   The  normalized 
velocity  profiles  of  the  pulsating  jet  downstream  from  the 
potential  core  conform  closely  to  the  probability  curve. 

This  test  was  conducted  in  the  fein  Engines  Labora- 
tory of  the  Mechanical  Engineering  Department  at  the  Univer- 
sity of  Minnesota,  Minneapolis,  Minnesota, in  partial  ful- 
fillment of  the  requirements  for  the  degree  of  Master  of 
Science. 
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Il\fTH0DUCTI01T 

Statement  of  the  Protlem: 

Many  studies,  "both  analytical  and  experimental, 
have  "been  made  concerning  the  mixing  in  fluid  streams.   The 
greater  majority  of  these  investigations  have  been  limited 
to  free  jets  and  steady  flow  conditions. 

Prandtl's  mathematical  analysis  (1)  of  mixing 
length  for  turbulent  flow  in  free  jets  as  applied  by  Tollmein 
(2)  along  with  Taylor's  theory  (3)  have  a  great  utility  in 
studying  various  types  of  fluid  jet  problems.   These  analyses 
do  not  offer  a  complete  explanation  of  the  problem. 

Squire  and  Trouncer  (4)  ha.ve  assumed  a  cosine  pro- 
file for  the  normalized  velocity  end  concentration  profiles 
dovmstream  of  the  potential  core.  For stall  and  Shapiro  (5) 
have  deduced  several  simple  empirical  formulae  for  the  mixing 
length  dovmstream  of  the  potential  core.  The  University  of 
Illinois  (6)  has  concluded  that  the  data,  of  free  turbulence 
can  be  well  correlated  by  probability  functions  and  applied 
Heichart's  Hypothesis  to  the  problem  of  mixing  length. 

Little  experimental  or  anal;/"tical  v;ork  has  been 
done  in  the  case  of  non-stationary  flov;  in  which  piilsations, 
either  rsjidom  or  regular,  exist.   The  fundajnental  equations 
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of  fluid  motion  may  "be  vn:itten  but  remain  unsolved  since 
pressure,  velocity  and  temperature  are  in  a.  continual  state 
of  rapid  random  fluctuations.   It  is  therefore  necessary 
that  experimental  analyses  be  made  which  might  introduce 
techniques  to  aid  in  analytical  treatment  of  the  problem. 

There  is  a  tremendous  need  for  design  data  con- 
cerning the  mixing  conditions  of  fuel  vapor  and  air.   Con- 
sideration has  been  given  to  turbulent  mixing  in  free  jets 
and  also  to  jets  issuing  into  ducts.   In  actual  conditions, 
such  as  pulse jets,  turbo jets,  ramjets,  and  combustion  cham- 
bers, pulsations,  both  regular  and  random,  do  occur.  Burton 
(8)  has  investigated  a  very  low  frequency  (250cpm)  pulsating 
jet.   In  this  investigation  close  agreement  was  found  of 
the  pulsating  flov/  with  the  stead;;/'  flow  empirical  formulae 
as  deduced  by  Forstall  ajid  Shapiro  (5). 

Godsey  and  Young  (7)  state  that  in  gas  turbine 
combustion  chambers  the  most  crucial  frequencies  are  in  the 
6000  cps,  250  to  600  cps,  and  25  to  60  cps  regions.  It  has 
also  been  found  that  the  critical  frequencies  in  the  pulse- 
jet  are  in  the  40  to  300  cps  region.  This  frequency  varies 
v;ith  the  size  of  the  pulsejet.  In  raj:njets  severe  pressure 
and  vibrational  fluctuations  are  encountered  because  of  the 
unsteady  nature  of  the  combustion. 
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Thus,  as  a  second  step  to  that  of  Forstall  8.nd 
Shapiro  (5)  it  is  necessary  to  introduce  pulsations  into 
the  floi\'  Tsefore  the  study  of  actual  conibustion  is  intro- 
duced. 

It  is  the  purpose  of  this  investigation  to  show 
the  extent  of  the  mixing  region  of  a  pulsating  air  jet  at 
varying  frequencies  in  a  steady  secondary  air  flow. 

Method  of  Approach: 

Basically  the  approach  to  this  investigation  is 
experimental  in  nature.   The  axial  velocity  field  is  plotted 
from  dynamic  pressure  obtained  hy  making  a  survey  of  the  mix- 
ing region  using  a  tota.1  head  tuhe  and  static  pressure  tap, 

A  method  of  measuring  varying  dynamic  pressures 
is  the  Statham  Low  Eange  Pressure  Transducer.   The  particu- 
lar model  availahle  for  this  experiment  had  a  natural  fre- 
quency of  750  cps.  The  Strain  Analyzer  and  Brush  Recorder 
units  record  dynamic  pressures  up  to  100  cps,  and  direction 
?s  well  as  magnitude  of  the  measured  pressure  can  "be  read 
from  the  chart. 

According  to  Prandtl  (1)  and  Tollmein  (2)  the  static 
pressure  is  constant  across  the  jet  within  small  limits. 
Therefore,  a  static  pressure  tap  was  located  in  the  side  of 


the  test  section  and  used  for  allL  points  within  the  flow  for 
each  particular  lateral  traverse. 

From  the  data  o'btained  a  comparison  of  the  mixing 
region  of  a  pulsating  jet  at  varying  frequencies,  including 
resonant  frequency,  msy  "be  made  with  the  mixing  region  of  a 
steady  jet. 
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TEST  EquiPI€ElJT 

General  Description: 

A   schematic  diagram  of  the  test  equipment  is  shovm 
in  Figure  1.  Figures  la,  ITd,  Ic,  Id,  le  are  photographs  of 
the  various  coraponent  parts  of  the  equipment  used. 

The  spparatus  consisted  of  two  separate  air  supply- 
systems  capa'Dle  of  furnishing  air  to  the  primary  and  secondexy 
sir  flov;  lines.   The  primary  air  flov/  was  ohtadned  from  the 
compressed  air  line  of  the  "building.   This  system  has  a  ca- 
pacity of  30  cul)ic  feet  at  a  pressure  of  100  p.s.i.  The 
primary  air  flov;  into  the  test  apparatus  was  manually  con- 
trolled "by  a  two  inch  gate  valve.  Prom  the  gate  valve  the 
air  passed  through  a  surge  tank  (2-^'  x  7'),  to  a  2"  stpjidard 
galvanized  pipe  which  contained  a  standard  A.S.M.E.  scua.re 
edged  orifice  v;ith  "radius"  taps  as  described  in  (10)  for  the 
measurement  of  the  mass  rate  of  floi'?  and  a  thermocouple  up- 
stream of  the  measuring  orifice  to  measure  temperature.   The 
primary  air  flow  v;as  then  led  into  ajiother  surge  tank  {2—*   x 
7')  and  then  into  a  2"  pipe  to  a  rotating  type  butterfly  valve 
v/hich  served  as  a  source  of  pulsations.  From  here  the  air 
was  smoothly  reduced  to  a  one  inch  inside  diameter  brass  tube 
and  ejected  into  the  test  section  as  a  free  jet.   The  cylindrical 
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surge  tanlcs  separated  the  butterfly  valve  froni  the  measuring 
orifice  to  prevent  pulsations  from  traveling  upstream  to  af- 
fect the  accuracy  of  the  measuring  orifice. 

The  source  of  secondary  air  was  a  gasoline  pov/ered 
Lycoming  Model  0-4o5-T  air  cooled  Army  tanlc  engine,  rated 
at  162  h.p.  at  2800  r.p.m.  which  drove  a  centrifugal  compres- 
sor.  The  compressor  was  a  7,48:  1  gear  ratio  supercharger 
taken  from  an  Allison  Y-1710  aircraft  engine.   The  speed  of 
the  hlov/er  was  controlled  by  throttling  the  tank  engine.   The 
air  from  the  compressor  vra,s  led  to  a  6"  I.D.  smooth  black 
iron  pipe,  through  a  measuring  orifice  as  described  in  (10) 
into  a  22"  x  35"  steel  drum  which  served  as  a  plenum  chamber. 
This  plenum  chamber  surrounded  the  primary  air  flow  line. 
From  the  plenvun  chamber  the  secondary  a-ir  passed  into  a  square 
bell  v.'ith  an  entrance  dimension  of  11"  x  11"  and  tapered 
smoothly  to  the  6  3/l6"  square  test  section.   Thus,  the  se- 
condary air  flov/ed  smoothly  into  the  test  section  s.nd  was 
concentric  to  the  primary  air  flov/.  An  18  mesh  wire  screen 
v/as  also  placed  14"  upstream  of  the  flov/  into  the  test  section 
to  attempt  to  get  sma-ll  sca,le  turbulence  liaving  the  same  pro- 
perties in  all  directions  at  the  test  section  entrance.  A 
thermocouple  wa-s  installed  in  the  seconda,ry  air  line  upstream 
of  the  orifice  to  measure  temperature  of  the  seconde.ry  air. 
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The  test  section  v;as  an  84"  long  square  (6  3/l6"  x 
6  3/16")  duct  made  of  20  gage  sheet  metal.   One  side  of  the 
duct  was  made  to  slide  in  sections  (14").   The  section  with 
the  total  head  tuhe  and  static  orifice  was  4"  wide.   Thus, 
this  section  coxild  "be  placed  in  any  position  along  the  test 
duct  longitudinally.   The  dovmstream  end  of  the  duct  was  open 
to  the  atmosphere.  A  nichrome  v?ire  was  centered  in  the  test 
section  such  that  a  circuit  was  completed  when  the  total  head 
tuhe  touched  it.   This  established  a  reference  point  for  cen- 
terline  values  and  the  start  of  lateral  traverses. 

The  total  head  tube  and  ststic  pressure  orifices 
were  connected  to  a  Stathajn  Low  Range  Pressure  Transducer, 
Model  P6-4D-2100  whose  natural  fre_quency  vras  750  c.p.s.   The 
gage  was  in  turn  connected  to  a  Strain  Analyzer,  Model  Bl-310, 
v/hich  led  to  one  side  of  a  double  pen  brush  recorder,  Model 
B 1-202. 

The  breaJcer  points  mounted  on  the  drive  shaft  of 
the  pulsator  were  led  to  a  brush  recorder  amplifier,  Model 
Bl-905,  which  in  turn  indicated  on  the  other  side  of  the  dou- 
ble pen  brush  recorder.  Model  Bl-202. 

Pulsator: 

A  2"  diameter  disc  type  butterfly  valve  rotating 
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in  the  primary  air  flow  pipe  provided  the  source  of  p\il  sat  ions, 
V/hen  the  valve  v;as  in  the  closed  position  the  2"  pipe  vras  com- 
pletely closed,  thus  the  air  flow  v;as  varied  from  zero  to 
majcimum.   The  "butterfly  valve  was  driven  by  a  direct  drive 
shaft.   This  shaft  had  a  standard  multiple  pulley  sheave  mounted 
v;hich  was  driven  through  a  system  of  tv/o  6"  variable  speed 
pulleys.   The  varia.hle  speed  pulley  system  vms  driven  through 
a  standard  multiple  pvilley  sheave  vrhich  was  attached  to  the 
shaft  of  a  ^r  hp,  110  volt,  1750  r.p.m. ,  A  C  electric  motor. 
Thus,  B   wide  range  of  frequencies  could  he  obtained  for  the 
pulsation  of  the  primary  air  flov;.  A  strobotac  was  used  to 
attain  and  accurately  maintain  the  desired  pvGLsating  frequency, 

Aajeccfiaatafcixj  jcam  v/as  moimted  on  the  pulsator  drive 
shaft  which  operated  a  set  of  breaker  points.  The  brealcer 
points  were  set  to  open  with  the  valve  in  the  fully  closed 
position..  In  this  manner  a  reference  point  for  valve  position 
was  obtained  on  the  binish  recorder  data. 

Instrumentation: 

A  Kiel  type  total  head  tube  vras  used  to  measure 
total  pressure.  The  venturi  shield  surroiinding  the  tube  tip 
to  insure  flow  normal  to  the  tube  opening  v;as  .375"  O.D.  , 
The  tube  opening  was  .0625"  O.D.  and  the  length  of  the  shrft 
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to  the  Statham  Lov/  Ra.nge  Pressure  Tra^nsducer  v;as  9  inches  long 
with  an  O.D.  of  .185".   The  Kiel  tube  was  mounted  in  the  4" 
sliding  section  of  the  test  duct  throu^  a  leak  proof  packing 
gland.   The  4"  sliding  section  also  contained  a  l/8  inch 
dipjneter  static  orifice. 

Vfater  filled  manometers  were  used  to  measure  the 
differential  pressures  across  the  square  edged  measuring  ori- 
fices in  the  2"  and  6"  flow  lines.  Upstream  static  pressures 
v/ere  measured  with  mercury  filled  manometers.   The  pressure 
taps  from  the  orifices  were  located  in  accorda.nce  v^ith  A.S.H.E, 
standards  for  "radius"  taps,  (10),  the  upstream  tap  being 
located  1  diameter  from  the  upstream  face  of  the  orifice  and 
the  downstream  tap  h   diameter  fron  the  dov/nstrea.ra  face  of  the 
orifice.   The  sta^tic  pressure  holes  were  :^-"  diameter,  free 
of  burrs  with  slightly  rounded  edges.   The  orifices  v/ere 
machined  to  an  I.D.  of  1.008"  and  4.002"  for  the  2"  prima.ry 
and  6"  seconda^ry  lines  respectively. 

Temperatures  of  the  orifices  were  measured  by  iron- 
constantan  thermocouples  inserted  in  the  flow  upstream  of  the 
orifices  in  accordance  v/ith  A.S.k.E.  standards.   Temperature 
readings  were  made  on  a  Brovra  potentiometer. 
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TEST  PROCEDUHE 

It  v/as  necessary  to  make  a  static  calibration  of 
the  Statham  Low  Range  Pressure  Transducer  using  the  Strain 
Analyzer  and  Brash  Recorder  instruments.  After  sllov;ing 
the  instruments  time  to  fully  v;arm  up  p  series  of  known  sta- 
tic pressures  were  applied  to  the  Statham  gage  and  the  deflec- 
tion read  on  the  Brush  Recorder.   The  scale  factor  on  the 
Strain  Analyzer  was  also  varied  from  200  to  50  micro  inches 
per  line  to  obta.in  the  desired  range  of  pressures.   Calibra- 
tion of  the  Statham  pressure  gage  is  shown  in  Fig.  2.   This 
calibration  was  used  to  reduce  all  data  obtained. 

All  the  instriimentation  was  assembled  and  checked 
for  lealcage  and  accuracy.   Prior  to  starting  e  ran  the  t3.nlc 
engine  was  allowed  to  warm  up  as  v;ere  the  Strain  Analyzer  and 
breaker  point  amplifier  unit.   Zero  points  v;ere  recorded  on 
the  Brush  Recorder  for  the  various  scale  settings. 

The  tank  engine  which  supplied  the  secondary  air 
flow,  was  brought  up  to  the  desired  speed  to  obtain  a  velo- 
city of  80  feet  per  second  in  the  test  section.   Then  the 
pulsstor  unit  v:as  turned  on  and  checked  with  the  Strobotac 
at  the  desired  frequency  setting.   iTov/  the  primary  air  flov; 
was  brought  to  the  desired  level  which  was  a  velocity  of 
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a.'bout  160  feet  per  second  at  the  exit  of  the  1"  "brass  tulje 
in  the  test  duct. 

Brush  recorder  data  v;as  then  taken  of  the  various 
traverses  of  the  total  head  tuhe.   Lateral  traverses  were 
made  at  a  series  of  stations  commencing  at  the  nozzle  exit 
and  continued  to  a  distance  of  72  nozzle  diameters  dov;nstrea.in. 
These  lateral  traverses  commenced  vrith  the  total  head  tuoe 
on  the  centerline  of  the  nozzle  and  continued  to  a  distance 
of  about  2.0"  from  the  centerline. 

This  type  of  r^on  was  made  vrith  the  pulsator  opera- 
ting at  three  frequencies:   1700  r.p.m. ,  1968  r.p.ra.,  and 
2200  r.p.ra.   The  frequency  of  1968  r.p.m.  vv'as  the  resonant 
frequency  of  the  system  as  found  experimentally.  A  plot  of 
the  resonant  frequency  determination  is  shown  in  Fig.  3. 

A  steady  flow  run  was  also  made  vrith  the  pulsator 
cut  off  and  the  "butterfly  valve  secured  in  the  full  open 
position.   This  was  made  for  comparison  v;ith  the  pulsating 
runs.   The  same  velocity  ratio  was  held  on  all  the  runs. 
Temperatures  end  pressures  v/ere  continually  checked  on  all 
runs. 


iS-Xt>VftlJ    x^'i 


13 


ANALYSIS  OF  BESULTS 

General: 

A  total  of  four  different  runs  v;ere  made  at  the 
same  velocity  ratio  as  follows:   (1)  Steady  State;   (2)  Re- 
sonant Frequency 1968  r.p.m. ;   (3)  2200  r.p.m.;  and  (4) 

1700  r.p.m.   The  dcta  obtained  is  listed  in  i'p'bles  I  to  IV 
inclusive. 

All  da.ta  obtained  from  the  recording  instrumenta- 
tion was  in  the  form  of  dynajnic  pressure  (ihs./in.^) .   This 
v/as  converted  to  velocity  (ft. /sec.)   as  shovm  in  the  section 
on  Sample  Calculations.   In  plotting  the  da,ta  the  nozzle  exit 
is  taken  as  the  origin  of  coordinates.   The  distance  along  the 
nozzle  centerline  is  denoted  hy  X/D,  positive  downstrea,m  ,  in 
nozzle  diameters,   Lrteral  distance  from  the  nozzle  center- 
line  is  denoted  "by  Y/D,   Thus,  the  station  at  the  nozzle  exit 
is  X/D=  0;  and  the  boundary  of  the  nozzle  is  Y/D  Z   .5. 

The  resonant  frequency  of  the  system  was  found  to  be 
1968  r.p.m.   It  was  necessary  to  find  this  experiment?,lly 
since  the  system  was  between  the  type  v;hich  has  one  end  open 
and  one  end  closed  and  the  type  v/hich  has  both  ends  open. 
Calculation  provided  a  starting  point  for  the  search  of  re- 
sonant frequency  determination.   Fig.  3  shov/s  a  plot  of  the 
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resonant  frequency  determination.   This  is  not  a  sharp 
resonejat  peak,  "but  one  of  a  damped  characteristic, 

i.  "butterfly  valve  provided  the  source  of  pulsp- 
tions  in  the  primary  air  flov;.   The  primary  flovr  varied  from 
zero  to  maximum  and  had  a  shape  similar  to  that  of  a  haT- 
monic  wave. 

If  flov;  reversals  occurred  the  total  head  tuhe 
v/ould  not  hfive  been  suitable  to  use.   At  no  time  did  the 
dynamic  pressure  approach  zero  on  the  recordinjf;  instrumen- 
tation, thus,  the  possibility  of  llovj   reversal  v/as  discounted. 

Steady  Flov/: 

It  was  desired  to  compare  the  mixing  region  charac- 
teristics of  the  tj'pe  of  flov:  actually  att^iined  with  those 
defined  by  previous  investi^p^ations.   Forstall  and  Shapiro  (5) 
have  presented  several  empirical  relations  for  isothermal 
jet  flows  similar  in  nature  to  this  experiment.   The  follov/ing 
axe  their  conclusions: 

(1)  "The  fully  normalized  shapes  of  the  velocity 
and  concentration  profiles  are  substantially  nlike. 

(2)  The  fully  normalized  shapes  of  the  profiles 
are  substantially  independant  of  velocity  ratio  and  of  axial 
distance. 
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(3)  Beyond  the  end  of  the  potential  core,  center- 
line  values  of  velocity  and  concentration  decrea.se  in  inverse 
proportion  to  increase  in  axial  distance." 

The  follov/ing  empirical  relations  are  "based  on  their 
experimental  data: 

(1)  X/D  value  for  the  end  of  the  potential  core: 

L  =  4  +  12 X 

(2)  Velocity  decay  downstream  of  the  potential 
core.  (X/D  >  L). 

U  -  U,  _  L 
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Fig.  4  is  p  plot  shov/ing  the  experimental  and  theo- 
reticpl  results.  Excellent  agreement  v;e.s  obtained  v;ith  the 
empirical  rate  of  velocity  decay.   This  is  indicated  "by  es- 
sentially ?  1:1  slope  of  the  logerithjuic  plots.   The  end  of 
the  potential  core,  "L" ,  is  found  hy  the  intersection  of  the 
stre.i£4"it  line  velocity  plots.   Using  the  empirical  formula 
an  L  =  9.85  is  obtained.   This  agrees  almost  exactly  v;ith 
that  plotted  experimentally. 

normalized  velocity  profiles  are  shovm  in  Fig.  9 

for  the  X/D  =  11  and  X/L  =  18  position.   This  was  done  by 

slotting  ^  ,  "  ^s  vs.  r/r^jj^.   The  values  of  r^^     are  ob- 

'^c  -  ^s 
tained  from  Fig.  12  v;hich  is  a  plot  of  Y/D  versus  velocity 

for  the  particular  X/D  position.   When  fully  normalized  in 


;Iqt  Bail 
:iwoIIo^   -•  . 


16 


this  manner,  all  profiles  downstream  of  the  potential  core 
are  almost  identical  for  the  steady  flov/  case.   These  pro- 
files have  a  close  resemblance  to  the  cosine  curve  as  found  "by 
Squire  and  Trouncer  (4).  For  direct  comparison  a  cosine  curve 

of  the  form  -R "  ^s  =  ^  (l  +  cosTTr)   is  shown.   Good  agree- 

c    s     d  ^^mv 

m.ent  is  seen  except  near  the  edge  of  the  jet  v/here  experimental 
accuracy  in  determining  U  -  Ug  is  not  good. 

Therefore,  v^ith  the  good  agreement  from  the  compari- 
sons v,'ith  previous  acceptahle  investigations  it  vas  thought 
that  the  accuracy  of  measurement  and  oversll  experimental  ac- 
curacy of  the  test  apparatus  was  sufficient  to  continue  with 
the  pulsating  flows. 

Pulsating  Flow: 

The  fiandsjnentaJ.  concern  of  this  investigation  v/as 
the  mixing  of  a  pulsating  air  jet  at  varying  frequencies  in  a 
steady  eir  flow.   Three  frequencies,  resonant  (1968  r.p.m.), 
one  ahove  resonance  (2200  r.p.m.),  and  one  belov/  resonance 
(1700  r.p.m.),  v/ere  used. 

Fig.  5  is  a  plot  of  the  centerline  values  of  the 
mean  velocity  for  the  pulsating  flov;  at  1968  r.p.m.  versus 
the  X/D  position.  It  is  seen  from  the  comparison  with  the 
empirical  ites.dy  state  formula  that  the  potential  core  is 
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shortened  considerably.   The  slope  of  the  velocity  decay- 
line  also  differs  from  the  steady  state, 

A  plot  of  the  centerline  values  of  the  mean  velo- 
city for  the  pulsating  flow  at  1700  r.p.m.  and  2200  r.p.m. 
versus  X/D  position  is  given  in  i'ig,  6.  Here  it  may  he  seen 
that  the  potential  core  is  shortened,  but  not  as  much  as 
that  in  the  resonant  frequency  case.   The  slopes  of  the  ve- 
locity decay  line  are  approximately  the  same  as  that  of  the 
resonant  frequency  run. 

The  following  is  a  table  to  enable  a  better  com- 
parison for  the  pulsating  runs  with  the  theoretical  steady^ 
state  ca,se: 


L  Slope 

Steady'  State  (Theor.)   9.89  ^/x/x, 

1968  r.p.m.            6.5  ^/(X/D)^-'^ 

1700  r.p.m.            8.8  ^/(x/D)^-^ 

2200  r.p.m.            8.0  ^/(x/D)^'^ 


For  the  steady  flow  care  velocity  decay  dovmstreejn 
of  the  potential  core  (X/D  >  L)  is 


^  "^s  -   L 
Up  -  U3    W 
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This  empirical  formula  for  the  pulsating  flow  for 
all  freouencies  appears  to  "be 


U  -Us  ~=   L   .  . 


From  the  above  comparisons  it  is  shown  that  pulsa- 
tions shorten  the  potential  core  and  give  faster  mixing. 
Resonent  frequency  gives  the  shortest  potential  core  and 
"best  mixing. 

Plots  v;ere  also  made  of  the  centerline  values  of 

the  amplitude  velocities  v/here  the  amplitude  velocity  is  the 

velocity  fluctuation  computed  from  the  meximum  and  minimum 

values  of  q,  versus  X/D  as  shown  in  Fig.  7  and  Fig.  8.   The 

BJnplitude  velocities  are  expressed  in  dimensionless  form 

0^  Uax  -  Us  •  From  these  plots  it  appears  that  the  cen- 

Up-  -Us 
terline  amplitudes  are  linear  with  viscous  damping.   The  slope 

of  the  resonant  frequency  run  is  less  than  thet  of  the  1700 

r.p.m.  and  2200  r.p.m.  runs.  Using  Fig.  7  the  equation 

U„      -.032  X/D   ^    , 
_  =  e  for  the  resonant  irequency. 


Ua  -U, 


The  plot.  Fig.  8,  for  17C0  r.p.m.  and  2200  r.p.m.  gives  a 
common  slope  for  "both  frequencies  and  the  equation 

Uax  -  ^s  ^  e-°^-^"  ^/^  . 

U     U 
^a  -  ^s 
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These  are  the  equations  for  the  decay  of  centerline  amplitude 
velocities  in  this  particular  test  apparatus.   Decay  of  the 
resonant  frequency  run  was  not  as  fast  as  that  of  the  1700 
r.p.m.  and  2200  r.p.m.  riins. 

An  attempt  vras  made  to  do  the  same  v;ith  the  lateral 
traverses  as  v/as  done  above  v;ith  the  centerline  amplitude 
velocities.   Hov;ever,  these  plots  were  non  linear  in  form. 
This  non  linearity  probably  is  caused  by  the  mixing  phenomena 
of  the  pulsating  primary  flow  and  the  steady  secondary  flow. 

The  data  shov/s  a  definite  phase  shift  for  center- 
line  values  of  X/D.   There  also  appears  to  be  a  phase  shift 
in  the  lateral  traverses  made.   However,  due  to  insufficient 
time  this  pliase  shift  v/as  not  recorded.   The  data  also  shov/s 
that  the  peak  velocity  through  the  butterfly  valve  did  not 
occur  st  the  full  open  position.  The  peak  velocity  occurred 
a,bout  10°  to  15°  from  the  fully  closed  position.   Apparently 
here  the  dyncjuic  pressure  was  at  its  highest  level. 

Again  the  profiles  v;ere  made  dimensionless  ^ni 

normalized  by  plotting  z. ~.  ^s  versus  r/r_,,  ,   The  re- 

suiting  profile  shapes  for  the  pulsating  frequencies  v/ere 
compered  v/ith  the  cosine  c\irve,  three-halves  pov;er  curve, 
and  error  curve.   The  curve  hnving  the  closest  resemblence 
of  the  above  three  v/as  that  of  the  error  curves.  According 
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to  He i chart ' s  Hypothesis  (6)  data  of  free  turbulence  can 

"be  correlated  by  probability  functions.   Thus,  since  the 

pulsating  flov/s  investigated  in  this  experiment  are  closely 

approximated  oy   the  error  curve  an  attempt  was  made  to  find 

the  curve  which  followed  the  trend  of  points  best.  A  curve 

of  the  form     U  -  Ug    [  l'^  (^     \2.4   v/as  found  to  be 

UT^  ""  [21    ^mv 
the  one  best  fitted.   This  curve  is  plotted  in  ?igs.  10  and 

11, and  again  it  is  observed  that  the  points  seem  to  spread 

near  the  base  of  the  profile.   The  experimenta.1  accura^cj'  in 

determining  (U  -  Ug)  at  the  edge  of  the  jet  is  not  too  good. 

The  regions  which  offer  the  larger  pressure  differentials 

give  better  results  than  those  near  the  edge  of  the  jet  v/here 

pressure  differentials  are  smaller. 

It  is  to  be  noted  that  the  r„,  value  reverses 

mv 

direction  between  different  X/D  stations  for  each  of  the  four 
runs  made.  The  r^^  value  is  defined  as  ths.t  radius  inhere  the 
velocity  is  the  aritiimetic  average  of  its  value  on  the  center- 
line  ajid  in  the  seconda.ry  stream  for  the  particular  X/D  sta- 
tion. For  the  ste8,dy  state  run  the  r^^^  value  reverses  be- 
tween X/d  =  24  and  X/D  ~   36.  For  the  pulsating  flows  the  rmv 
value  reverses  direction  as  follov/s: 

Resonant  frequency  (1S68  r.p.m.)  — X/D=  15  and 

X/D  =  18 
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1700  r.p.m. X/D  -   12  and  X/D  =  15 

2200  r.p.m. X/D  =  18  and  X/D  »  24 

Twice  the  Tjjj^  value  defines  the  radius  of  the 
outer  edge  of  the  velocity  profile.  The  outer  edge  of  the 
velocity  profile  never  reaches  the  sides  of  the  test  duct. 

For  the  pulsating  flows  the  rjjj^  value  reverses 
direction  further  downstream  as  the  pulsating  frequency  in- 
creases.  The  cause  for  this  was  not  ascertained  although 
it  is  thought  that  phase  shift  might  he  a  contributing  factor. 

Points  at  which  the  centerline  flov;  might  "be  con- 
sidered mixed  to  various  degrees  are  shovm  in  the  follov/ing 
ta'Dle.   Tliese  values  are  obtained  from  Figs.  4,  5,  and  6. 
The  percent  mixing  is  based  on  1 1  -fu  -  Ug\>.100,  v/here 


I   LUp  -  U^) 


U  -  Ug  is  the  sctual  velocity  gradient  v/hich  exists  at  a 
particular  point  and  U^  -  Ug  is  the  maximum  velocity  grg.- 
dient  which  can  exist.   Thus,  the  mixing  is  complete  when 
the  fraction  U  -  Ug  reduces  to  zero. 


Up  -  Ug 

steady  Flow 

X/D  = 

50)6 
19 

75'^ 
37 

90% 
90 

He  sonant  Freauency 
(1968  r.p.m.) 

X/D  = 

10.8 

.19 

42 

1700  r.p.m. 

X/D  = 

15.5 

29 

66 

2200  r.p.m. 

X/D  = 

13 

24 

54 

22 


From  the  above  table  it  is  shown  that  pulsating 
flov;  mixes  more  rapidly  than  stesdy  flow.   It  is  also  seen 
that  at  resonant  frequency  the  mixing  is  accomplished  in 
about  one  half  the  distance  recuired  by  steady  flow. 
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COHCLUSIOIIS  MD  RECOMIffilTOATIONS 

1.  The  velocity  mixing  region  of  an  isotherm?! 
pulsating  jet  is  much  shorter  thaji  tha.t  of  a,  steady  flow  jet 
of  the  S8.iiie  configuration  at  the  same  velocity  ratio. 

2.  Pulsations  at  the  resonant  frequency  of  the 
test  apparatus  gives  the  shortest  potential  core  and  the 
"best  mixing  chs.racteri sties.  At  frequencies  above  and  helov/ 
resonance  the  potential  core  is  shorter  than  the  stea.dy  flow 
case  and  the  mixing  is  better. 

3.  The  centerline  velocity,  dov.mstreain  of  the 
potential  core,  decays  in  direct  proportion  to(X/D)  *   for 
the  pulsating  flovrs. 

4.  The  decay  of  the  centerline  amplitude  veloci- 
ties appears  linear  with  viscous  da.mping.   The  decay  of  the 
resonant  frequency  run  is  not  as  fast  as  that  of  the  runs 
above  and  belov/  resonance,   ?or  the  resonant  run 

^ax  -  Ug  _   -.032X/D    ig  the  ecuation  of  the  decay. 


34  - 


The  eouation  for  the  decay  above  and  helow  resonance  is 

■ re*     '  ,  Decay  of  the  lateral  traverses  for 

all  pulsating  runs  was  non-linear. 


5,  The  normalized  velocity  profiles  dov/nstream  of 
the  potential  core  are  essentially  of  the  same  shs,pe,  irre- 
spective of  the  X/D  position.   They  closely  resemble  the 
probability  curve  more  so  then  the  cosine  or  three  halves 
pov/er  curve.  A  curve  of  the  form   U  -  Us   flu^  1  '"* " 

fits  the  points  v^ell. 

6,  Good  agreement  was  obtained  vdth  the  findings 
of  Forstall  and  Shapiro  for  the  steady  flov/  case.   The  loca- 
tion of  the  potential  core  is  closely  defined  by  the  empirical 
relation  L  =  4  +  12  A.   The  centerline  velocity,  downstream 

of  the  potential  core,  decays  in  direct  proportion  to  X/D. 
The  normalized  velocity  profiles  closely  resemble  a  cosine 
curve. 

7,  The  centerline  flow  was  considered  90>d  mixed 
at  X/D  positions  of  90  for  steady  flow,  42  for  resonant 
frequency,  66  for  1700  r.p.m. ,  and  54  for  2200  r.p.m.   There 
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was  75^t;  mixing  at  X/D  positions  of  37  for  steady  flov;,  19 
for  resonant  frequency,  29  for  1700  r.p.m,,  and  24  for  2200 
r.p.m, 

8.  With  respect  to  the  test  apparatus,  it  is  re- 
commended that  a  surge  tank  "be  placed  in  the  secondary  air 
line.   Small  amplitude  pulsations  were  noted  on  the  record- 
ing data  which  v/ere  traced  to  the  "blower  operated  by  the 
tank  engine. 

0,  Further  investigations  of  this  nature  are  re- 
commendeci  to  find  the  effects  of  phase  shift,  pulsation  form 
and  velocity  ratio  upon  the  mixing  region  of  s   pulsating  jet. 
Phase  shift  is  apparent  and  readable  to  a  low  degree  of  ac- 
curacy on  the  data  obtained.   It  is  thought  that  photographs 
of  an  oscilloscope  screen  might  provide  reada.ble  data. 
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FORl-OTJLS  MTD  SAI^IFLE  CALCULATIONS 


1,      Air  Flov;  Measurement  and  Calculation 

The  equation  used  to  find  mass  flow  wasi 

1 


V  =  0.688  As  Key  (f  ^-13)2 
(see  Ref.    10)     '      ' 

2,  Calculation  of  Reynolds  Kuciber 

%  =     P  u  D 

3,  Frequency  Computations 

The  natural  frequency  of  the  tube  v.dth  one  end 
closed  is  computed  from 

f  =  a 

4L 

The  natural  frequency  of  the  tube  with  both 

ends  open  is  computed  from 

f  =  a 
2L 

4,  Calculation  of  velocity 

p  =    p  RT 

where  p  =  2116  Ibs./ft^ 

^  =  1715  ft^/sec^  OR 
T  s   565^  R 

f21l6  Ibs./ft^  nnoiD    ^i-.,..^/^+3 

=      (1715  ft'^/sec^°R)(565->R)  =      '^^^^^   ''^^'/^' 
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q  =  -^V2 

where  f  =   .00218  slugs/ft^ 

q  =  Dynamic  Pressure  I'bs./in,^ 

V^  =  2a  =  (2) (144  in^/fts)(  q  Ibs./inS) 
"p  .00218  slugs/ft^ 

V^  =   132  X   10^  q  ft^ 
sec^ 

V  =  363.5  VT*   ft/sec. 
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NOffiHCLATUEfl 

D  diameter  of  primary  jet  noszle,  inches. 

L  value  of  X/D  for  end  of  potential  core. 

r  radial  distance  from  jet  axis. 

q  dynanic  pressure,  psi, 

Tj^^. . . .  radius  where  velocity  is  arithmetic  average  of  its 
value  on  the  axis  and  in  secondary  stream,  for  a 
given  X. 

U time  average  velocity  in  axial  direction, 

Uq velocity  on  axis  at  a  given  X, 

U  ....  velocity  of  primary  stresjn  at  X  =  0. 

Ug  ....  velocity  of  undisturbed  secondary  stream. 

"U^   ....  velocity  of  pulsation  amplitudes  on  the  axis  at  X  =  0. 

Ug,-  ...  velocity  of  pulsation  amplitudes  on  the  axis  at  a. 
given  X, 

X  axial  distance  downstream  from  nozzle  exit,  inches, 

Y  lateral  distance  from  noxzle  centerline,  inches. 

J^ velocity  ratio.  Ug/Up. 

p  density 

^  viscosity 
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APPENDIX 


X 

D 

Scale 
Factor 

0 

100 

3 

100 

6 

100 

7 

100 

8 

100 

9 

100 

10 

100 

11 

100 

12 

100 

15 

100 

18 

50 

24 

50 

36 

50 

48 

50 

60 

50 

72 

50 
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TABLE  I 

STEADY  FLOW  CEUTSRLIIEE  VELOCITIES 
Up  =  160ft/sec;  Ug  =  78  ft/sec;  A=  .487 

Deflection 
(m.m. ) 

8.8 

9.2 

8.7 

9.0 

8.5 

8.6 

7.3 

6.9 

6.3 

5.2 

8.3 

7.0 

6.3 

5.3 

5.4 

5.4 


bs/in.2 

U2  X  10^ 
ft^/sec^ 

U 
ft/sec 

U  -  Us 

.195 

25.2 

159.1 

.205 

26.4 

162.0 

1.025 

.194 

25.0 

158.1 

.98 

.200 

25.8 

160.5 

1.005 

.190 

24.5 

156.6 

.96 

.191 

24.6 

157.0 

.964 

.162 

20.9 

144.5 

.935 

.154 

19.88 

141.0 

.89 

.140 

18.10 

134.8 

.694 

.115 

14.82 

121.9 

.535 

.107 

13.8 

117.5 

.482 

.091 

11.72 

108.1 

.367 

.082 

10.59 

103.0 

.305 

,069 

8.90 

94.5 

.201 

.070 

9.03 

95.1 

.208 

.070 

9.03 

95.1 

.208 

-   33  - 

TA£LE  II 

STEADY  FLOV/  LATERAL  TRAVERSE  FOR  X/D  =  11 

Up  =   160  ft/sec;   Uc  =  141  ft/sec;    Ug  =  78  ft/sec;    A  =   .487 

r„^  =  .42" 
mv 


y 

D 

Scale 
Factor 

Deflection 
(m.m. ) 

lts/in2 

U^  X  10^ 
ft^/sec^ 

U 
ft/sec 

r 

U  -  Us 

0 

100 

6.9 

.154 

19.88 

141.0 

0 

1 

.1 

100 

6.7 

.150 

19.35 

139.0 

.238 

.959 

.2 

100 

6.3 

.140 

18.05 

134.5 

.475 

.850 

.3 

100 

5.6 

.125 

16.11 

127.0 

.714 

.70 

.4 

50 

8.7 

.113 

14.60 

121.0 

.955 

.57 

.5 

50 

7.8 

.100 

12.90 

113.8 

1.19 

.415 

.6 

50 

7.3 

.095 

12.25 

110.8 

1.43 

.35 

.8 

50 

5.8 

.075 

9.69 

98.1 

1.90 

.077 

1.0 

50 

5.4 

.071 

9.15 

96.6 

1.2 

50 

5.3 

.069 

8.90 

94.5 

1.4 

50 

5.3 

.069 

8.90 

94.5 

1.6 

50 

5.2 

.068 

8.77 

93.5 
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TABLE  III 

FULSATIEG  FLOW  CEIJIERLIIIB  MEAIT  AIJD  AMPLITUDE  VELOCITIES 

EESONAM!  PHEqUElICY  -   1968  r.p.m. 

Up  =   158  ft/sec;  Us  =  77.5  ft/sec;    A  =   .^9 


^ 

Factor 

(m.m.) 

(m.m.) 

lb/in2 

It/in^ 

ft/sec 

ft /sec 

Up  -Us 

Ua-Us 

0 

200 

4.5 

11.0 

.195 

.477 

160 

251 

— 

-.. 

3 

200 

4.1 

9.5 

.180 

.412 

154.5 

255 

.959 

.91 

6 

200 

4.0 

9.0 

.175 

.390 

152 

227 

.927 

.865 

7 

200 

3.5 

8.5 

.153 

.570 

142.1 

221 

.804 

.825 

8 

200 

3.2 

8.0 

.140 

.348 

136 

214.5 

.727 

.79 

9 

200 

2.8 

7.8 

.124 

.339 

128 

212 

.628 

.775 

10 

200 

2.5 

7.0 

.108 

.304 

119.5 

200 

.521 

.706 

11 

200 

2.3 

7.0 

.100 

.304 

115 

200 

.465 

.706 

12 

200 

2.15 

6.3 

.093 

.274 

111 

190.5 

.416 

.651 

15 

200 

1.85 

6.0 

.080 

.260 

103 

185 

.316 

.62 

18 

200 

1.70 

4.4 

.075 

.191 

99.5 

159 

.274 

.47 

24 

200 

1.45 

3.8 

.064 

.166 

92 

148 

.180 

.405 

36 

100 

2.6 

3.1 

.057 

.070 

87 

96 

.118 

.107 

48 

100 

2.4 

.055 

85.2 

.096 

60 

50 

4.0 

.051 

82.1 

.057 

72 

50 

4.0 

.051 

82.1 

.057 
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TABLE   IV 

FUISATING  FLOW  LATERAL  TEAVISHSE  FOR  X/D  =   11 

HESOKAHT  FEIDCiUEilCY  -   1968  r.p.m. 

Up  =  158  ft/sec;  Ug  =  115  ft/sec;   Ug  =  77.5  ft/sec 

X  -   .49;   r^„  =   .48 


y 

Scale 
Factor 

Deflec- 
tion 
(m.m.) 

lbs/in^ 

ft/ sec 

r 

Uc  -  Ug 

0 

200 

2.3 

.100 

115 

0 

1 

.1 

200 

2.25 

.097 

113 

.208 

.947 

.2 

100 

4.15 

.092 

110 

.416 

.868 

.3 

100 

3.9 

.036 

106.5 

.625 

.774 

.4 

100 

3.6 

.080 

102.6 

.833 

.670 

.5 

100 

3.0 

.066 

93.5 

.960 

.426 

.6 

100 

2.5 

.055 

85.1 

1.25 

.202 

.8 

100 

2.1 

.047 

79.0 

1.66 

.04 

1.0 

100 

2.0 

.045 

77.2 

2.08 

0 

1.2 

100 

2.1 

.047 

79.0 

1.4 

100 

2.0 

.045 

77.2 

1.6 

100 

2.0 

.045 

77.2 

1.8 

100 

2.0 

.045 

77.2 

>3\n  8c X  =  qU 
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15 

K 

^ 

^ 

i/> 
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v 

i 

\b 

vx 

? 

? 

3 

Hj 

« 

/ 


^ 


pO  (^ 


1     Bneine, 

it   set-up  SiiOwxi  is  uot  part  oi   t>^is  experiment) 


Fig,   l-T)     Test  duct  and  recording  instrumentation. 
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/          ^^^^  ^^            Infl            IL  \ 

Pig,   1-c     Palsator  Unit 


Fig,   1-d     Controlling  inso-ixwicaoaoxua. 


Fig.  1-e  Control  panel  for  en^ne  operation. 
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